Scientific Question: Given the significant reduction for the risk of breast cancer after obese patients undergo bariatric surgery, what specific genes are responsible for the reduced risk and how do the expression levels of these certain genes change after the surgery?
Hypothesis: If the levels of Estrogen Receptor Alpha decrease after bariatric surgery, then by analyzing the specific gene (ESR1) expression levels of this protein and associated pathways will show a marked change from pre-weight loss surgery levels that are responsible for the preventative phenotype against breast cancer.
The packages Edger and limma are used to filter out lowly expressed counts and to model the voom plot in RNA-seq and also EdgeR is used in Pathway analysis to gain access to cerrtain functions
library(edgeR)
## Loading required package: limma
library(limma)
The package msa used to align multiple sequences from a single fasta file with its function msa
library(msa)
## Loading required package: Biostrings
## Loading required package: BiocGenerics
##
## Attaching package: 'BiocGenerics'
## The following object is masked from 'package:limma':
##
## plotMA
## The following objects are masked from 'package:stats':
##
## IQR, mad, sd, var, xtabs
## The following objects are masked from 'package:base':
##
## anyDuplicated, append, as.data.frame, basename, cbind, colnames,
## dirname, do.call, duplicated, eval, evalq, Filter, Find, get, grep,
## grepl, intersect, is.unsorted, lapply, Map, mapply, match, mget,
## order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank,
## rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply,
## union, unique, unsplit, which.max, which.min
## Loading required package: S4Vectors
## Loading required package: stats4
##
## Attaching package: 'S4Vectors'
## The following objects are masked from 'package:base':
##
## expand.grid, I, unname
## Loading required package: IRanges
## Loading required package: XVector
## Loading required package: GenomeInfoDb
##
## Attaching package: 'Biostrings'
## The following object is masked from 'package:base':
##
## strsplit
seqdata<- read.csv("/Users/heatherotto/Downloads/GSE126803_AI007_Readcounts.csv")
seqdata